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1 Introduction 



Deep inelastic lepton - nucleoli scattering (DIS) has been an important tool in particle 
physics for more than twenty-five years. Recent high precision measurements by groups 
such as the NMC, SMC and various SLAC groups [|I], ^ |], |^ are testing the limits of our 
theoretical knowledge of the structure of the nucleon. The SMC and the E143 groups H 
have recently reported the first measurements of the transverse spin structure function 
G2{x, Q"^) [] , which has a twist-three component qt- In the near future we can also expect 
the HERMES experiment to investigate the spin dependent structure functions in some 
detail. Also the development of polarized beams at RHIC will lead to measurements of 
the chiral-odd distributions hi^i, which are twist-two and twist-three respectively 0. 

Higher twist structure functions will also be important at the energies of CEBAF and 
the proposed ELFE accelerator because of the phenomenon of hadron-parton duality § , 
or 'precocious scaling' |]10|]. This predicts that at even moderate values of Q^, higher twist 
contributions to observables, such as the nucleon form factor, will be non-negligible. Thus 
it is important to have some theoretical understanding on the expected sizes and shapes 
of the higher twist structure functions. 



In a recent paper, Ji |Tl| defined the 14 possible nucleon structure functions in the 
standard model, and calculated each structure function in terms of all the possible parton 
distribution functions up to the level of twist-four (ignoring loop effects from QCD radia- 
tive corrections). The parton distribution (or correlation) functions are defined in terms 
of Fourier transforms of matrix elements of non-local quark and gauge fields separated 



^Throughout this work I use capital letters F, G to denote the structure functions measured in exper- 
iment and lower case letters /, g, h to denote the various parton distributions. 
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along the light-cone [IT], |T^. At present these matrix elements cannot be calculated 
in QCD, however calculations of twist-two matrix elements using the wavefunction of the 
MIT bag model have been performed [0, and these describe current data relatively 
well. In this paper these calculations are extended to all the parton distribution functions 
of twist-three and four containing only two quark fields. 

2 Parton Distribution Functions 



Following Ji |TT[], let us consider a nucleon of mass M with momentum and polarisation 
vector S'^. If we choose the nucleon to be moving along the 2;- axis with momentum P we 
have 

Pf" = {VM^ + P\0,0,P). (1) 
We introduce two orthogonal, light-like vectors p and n 



pi^ = i(VM2-fP2 + p)(i^ 0,0,1), 

= -^(VM2 + P2_p)(i^ 0,0,-1), (2) 

satisfying = + We also decompose the polarisation vector, S = S\\ +MS±, 

where 

S^^=lf-^n^, 5^: = (0,1,0,0). (3) 

A parton distribution with k light-cone momentum fractions M{xi, . . . , x^) is defined 
via the matrix element 

/ n eM^\x,){PS\Q{X,n, . . . , Xkn)\PS) = M(xi, . . . , Xk)f{p, n, S^) (4) 



where Q is a product of k quark and gluon fields, and T is a Lorentz tensor. 

If the mass dimension of the parton distribution M(xj) is (Im and that of T is (It, then 
at large the matrix element will behave as A'^^Q'^^, where A is a soft QCD mass scale 
related to non-perturbative physics. A distribution with mass dimension (Im is called 
a twist- (c/a/ + 2) distribution. As the dimension of a physical observable is fixed, the 
higher (Im becomes, the higher the inverse power of hard momenta must become in the 
observable. Thus a twist-n parton distribution can only contribute to physical observables 
of twist-n or higher, which behave as Q^^'-"'^"^ where a > 0, in the Q"^ ^ oo limit. In 
the following we will often rescale a parton distribution by factors of (M/A) to make it 
dimensionless, but the twist wilU be unchanged. 

For scattering processes involving two quark fields we can define a quark density matrix 

/■ r/A 

M„^(x) = J —e'^^{PS\MO)M>^n)\PS). (5) 

From this it is possible to systematically generate the possible distributions at a given 
twist. The twist-2 part of the density matrix can be written 

^(^)ltwist-2 = \ ^/iW + ^75 ■ n)gi{x) + ^75 ^± ^hi{x) (6) 
where the three quark distribution functions /i, gi and hi are defined 

hi{x) = \j'^e^^^{PS^\mil^^±^{>^n)\PS^). (7) 

These represent the unpolarized quark density, the quark helicity density and the quark 
transversity density [0 respectively. The light-cone gauge, = 0, has been chosen, so 
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that the distributions are manifestly gauge invariant. 

Similarly at the twist-3 level, the appropriate portion of the quark density matrix can 
be written 

^(^)ltwist-3 = ^[e(^) + ^(-^ll • n)(ri fl- fi ^HhL(x) + 75 ^±9t{x)] (8) 
where A is a soft mass scale in QCD. The three twist-3 distribution functions are given 

by 

9t{x) = ^J^e^'^{PSM0H^±i^{Xn)\PS^). (9) 
At twist-4 we have 

A2 

^(^)ltwist-4 = -^W^i^)+ 1^^5{S\\ ■ n)g3{x)+ fl-f^ ^±h3{x)] (10) 
with the twist-4 distributions 

9s{x) = ^/^e^^^(P5|||V;(0)75#(An)|P5||), 

hsix) = 2^ / ^e^'^iPS^OH ^± m^n)\PS^). (11) 

The quark field -0 can be decomposed into 'good' and 'bad' components, ■0+ and -0- 
respectively 

= pH, + 7^ ^ 7!^^° ^ 

By inspection it can be seen that the twist-2 quark distributions involve only the 'good' 
component ■0+, whereas the twist-3 distributions involve mixing one 'good' and one 'bad' 
component, and the twist-4 distributions involve only the 'bad' components. 
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The QCD equations of motion |T5 



2^^_(An) = 1 ^(-zp^ + m,)V'+(An), (13) 
-i^^.{\n) = ]^^+{\n){-ip^ + m^) (14) 

make it possible to eliminate the 'bad' components from the twist-three and four distri- 
butions, at the cost of introducing gluon fields into the matrix elements. Because the 
model wavefunctions do not include gluon fields, we will not do this here. Also note that 
at twist-three and twist-four there exist distributions involving the gluon field with two 



or three light-cone momentum fractions, such as ||TT|, |12 



E{x,y) = f ^^e^^V^(^--)(P|^(0)^/)x(H ^^(An)|P), (15) 

and 

Biix,y,z) = ^g^e^^V^(^--)e^'^(^^^)(P|^(0) /)^(z.n)z ^^(H^(An)|P), 

(16) 

which can be related to e{x) and /^{x) by integrating over y or y and z respectively. 
However, as the MIT bag wavefunction has no explicit gluon field, these distributions will 
be zero in the model. Also the distributions with only one light-cone momentum fraction 
are of the most interest for DIS and Drell-Yan processes. 

Finally there exist distributions at the twist-four level involving four quark operators 

eg 

Uax,y,z) = -]- [ ^^^e^Ax.^^^(^-)e^'^(^-^)(P|^(0) ^^(z/n)^(/in)V;(An)|P), (17) 

which is a four quark light-cone correlation function. Calculating this distribution by the 
method below would require the evaluation of the overlap integral between the four quark 



fields over the bag volume, which is expected to much smaller than the two quark overlap 
integral required for the two quark correlation functions. Hence these will not considered 
further in this work. 

3 Calculation of Quark Distributions 

At present no QCD wavefunction for the nucleon can be calculated. So in order to make 
useful calculations of the quark distributions at any twist it is necessary to the use the 



wavefunction from some phenomenological model of the nucleon. The MIT bag model |T6 
is used here as it incorporates relativistic, light quarks, and also models confinement. It 
also has the further advantage that the wavefunction is simple and analytic. Other models 



could also be chosen 17 



The major problem in calculating the relevant matrix elements for the quark distri- 
butions is ensuring that momentum conservation is obeyed throughout the calculation, 
hence ensuring that the calculated distributions have the correct support, i.e. they vanish 
for light-cone momentum fraction x outside the interval [0,1] |T^. To guarantee mo- 
mentum conservation, a complete set of intermediate states, J2m \m){m\, can be inserted 
into the matrix elements of the quark distributions (eqns. (0, |lll)). Using translational 



invariance of the matrix element, all the n dependence can be moved into the argument 
of the exponential function. Then integrating over A gives a momentum conserving delta 
function. The twist-two quark distributions then become 



/i(^) 



-=E<5(p+(l-x)-p+)|(m|^+(0)|P)p 

V m 



g,ix) = i=E'5(Pni-^)-P^)[K^|i?^+(0)|P5||)P-|(m|LV^+(0)|P5||)p] 

m 



m 

|(m|g_V'+(0)|P5^)|^], (18) 

where R{L) is the projection operator for right-(left-) handed quarks R{L) = (1 ± 75)/2, 
and Q± is the projection operator Q± = (1 ± 75 ^±)/2, which projects out eigenstates 
of the transversely projected Pauh-Lubanski operator ^±^5 in a transversely projected 
nucleon. 

The twist-four distributions are similar to those of twist-two, except they involve 'bad' 
components of the quark wavefunction 
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f4{x) 



-=;E%+(l-a;)-p+)|(m|^_(0)|P)r, 
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^3(:^) = ^E%^(i-^)-p;:;)[iHi^V'-(o)|p^ii)r-i(m|it:v'-(o)|P5ii)r], 

h{x) = i=E%ni-^)-P+)[IHQ+V'-(o)|P^±)r- 

\{m\Q.i;.{Q)\PS^)\']. (19) 



Here the distributions have been rescaled to make them dimensionless, and factors of 
{M/Ay have been absorbed in the twist-four part of the density matrix, ^(a^)ltwist-4' 

Note that the twist-two and twist-four distributions have a natural interpretation in 
the parton model, where they are related to the probability of finding a parton carrying 
fraction x of the plus component of momentum of the nucleon, and in the appropriate 
helicity or transversity eigenstates. In the twist-three case, the distributions do not have 
a similar interpretation in the parton model. However we can still guarantee momen- 
tum conservation by introducing a complete set of intermediate states, J2m |^)("^|) 
then write the distributions in terms of the matrix elements between nucleon states and 



intermediate states 

e(x) = 5:5(p+(l-x)-p+)(P|^t(o)|^)^(/)^,(^|^(0)|P),, 

m 
m 

9t{x) = Y.^iP^i^-^)-Pt^){PS^\^PHO)\m)a{l'l' ^±)at{m\m\PS±)b (20) 

m 

where the matrix indices have been shown exphcitly. Again the distributions have been 
rescaled so that they are dimensionless, absorbing factors of M/ A into the density matrix 

^(^)ltwist-3- 

The next step in the model calculation of the distributions is to form the momentum 
eigenstates \P) and \m) from the static states of the model. This can be done using either 
the Peierls-Yoccoz projection, which gives a momentum dependent normalisation, or 
the Peierls-Thouless projection, which leads to a more difficult calculation, but which 



preserves Galilean invariance of the matrix elements. Using either method, the correct 
support for the distributions is guaranteed by the above formalism. The distributions 
can then be calculated in terms of the Hill- Wheeler overlap integrals between the quark 
wavef unctions |]TB|, |T3 . 



Using the wavefunction of a model also introduces a scale into the calculated distri- 
bution functions. This is the scale at which the model wavefunction is considered a good 
approximation to the true QCD wavefunction, which is presently unknown. The natural 
scale for the bag model, and most other phenomenological models employing light rela- 
tivistic quarks, is the typical transverse momentum of the quarks, ~ 400MeV. In order 
to compare a calculated distribution function with experiment, the calculated distribution 
needs to be evolved from the model scale up to the experimental scale Q^. This has previ- 



ously been done using leading order QCD evolution for the twist-2 distributions and 



gi{x) |jT3|, [T^, with good agreement being obtained for a value of in the region of 250 



500 MeV. This could be criticised on the grounds that the strong coupling constant is not 



small in this region, however calculations using next to leading order evolution I^T] also 
give good agreement with experiment for values of ~ 350MeV and as{Q'^ = /U^) ~ 0.75. 

For the higher twist distributions there do not yet exist comprehensive calculations 
of the relevant anomalous dimensions for evolution of the distribution functions. For the 
twist-three distributions qt and hi full calculations of the anomalous dimensions have 
appeared in the literature ^ For the lowest moments of these distributions 
[n = 3,4), the leading order evolution gives similar results to naive power counting for 
\/Q^ in the region of 1 GeV. Thus we shall neglect QCD corrections to the evolution of 
the higher twist distributions, and use naive power counting to evolve these distributions. 

In Figure 1 the nine twist-2, 3 and 4 distributions involving two quark correlation 
functions, calculated at the bag scale fi, for a bag radius of 0.8 fm are shown. The 
Peierls-Yoccoz projection has been used to form the momentum eigenstates \p) and |m), 
and only intermediate states containing two quarks have been considered. As is expected 
from equations (|18| , [l^, ^OD , each of these distribution functions are similar in magnitude 



to one another. It is worth noting that the parton distributions satisfy the equalities 

fi{x)+gi{x) = 2hi{x) 
e{x) + = 2gT{x) 

Mx)+g3{x) = 2h{x) (21) 

which are the lower bounds of Soffer's inequalities [p5| ]. 
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The evolution of the twist-3 distributions is particularly interesting as the distributions 
have contributions from local operators of both twist-3 and twist-2. Hence the moments of 
the twist-3 distributions are related to the moments of the twist-2 distributions 0, ^ |24 



Mn[e] = ^Mnlfl] + Mn[e^'^] 

^-^St] = ^Mn[9i]-l^Mn-i[h,]+Mn[9?] (22) 

where m is the quark mass, A^„[/] = Jq x^f{x)dx is the nth moment of the distribution 
f{x) and f^^^ denotes the genuine twist-3 part of the distribution f{x). Note that when 
contributions from operators proportional to the quark mass or from twist-3 operators 
are ignored, the third of these relations is equivalent to the Wandzura-Wilczek relation 



g,(x) = -g,{x)+ f'^dy. (23) 
Jx y 

In what follows we shall ignore terms directly proportional to the quark mass, as the 
masses of the u and d quarks are zero (or very small) in the MIT bag model. (If this 
discussion were extended to include strange quarks and anti-quarks then it would be 
important to keep terms in the quark mass.) 

We can invert the moments above to find the genuine twist-3 parts of each distribution 
in terms of the calculated distributions and the twist-2 distributions 

e(3)(x) = e{x) 

hf{x) = hL{x)-2x C%h,{y) 

Jx y^ 

gf{x) = grix)- f %iiy). (24) 
Jx y 
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This separation makes it possible to evolve the calculated distributions e,hL,gT from the 
model scale up to experimental scales. The procedure I use for the evolution is to separate 
the calculated distributions at the model scale = ji^ into twist-2 and twist-3 parts using 
eqn. (p^. The twist-2 parts are then evolved using the Gribov-Lipatov-Altarelli-Parisi 



(GLAP) equation to leading order, whereas the twist-3 parts are evolved according 
to naive power counting f^^\Q'^) ~ The model scale has been chosen as /i = 0.4 

GeV, and in Figure 2 the full twist-3 distributions at = 1 and 10 GeV^ are shown. 
Comparing between e{x,Q^), which has no twist-2 part, and hL{x,Q^) and gT{x,Q^) 
enables us to see how the twist-2 parts of the latter two distributions dominate at higher 

In fact it is probably not a particularly good approximation to evolve e{x, Q"^) accord- 
ing to naive power counting. While a calculation of the anomalous dimensions for the 
operators contributing to e{x) has yet to be done, we can surmise a few facts about the 
lowest moments, including that the lowest moment of e{x) will have a small anomalous 
dimension despite the corresponding operator being twist-3. The relevant local operators 
for e{x) in the operator product expansion are 

0WM2.../.n _ . . . {zD^'")iIj - traces (25) 

where iS„ symmetrizes over the Lorentz indices /xi, . . The relevant matrix elements 
of these operators are 

(^p\Qm-t^n^p^ = 2Me„P^ip^2 . . . P'^" - traces. (26) 

Standard techniques give the sum rules 

en = J dxx'^e{x). (27) 

12 



In particular, for n = we have 

/ dxe{x) = ^{P\^^l;\P) (28) 

which is related to the nucleon a term, (P||(mu + md){uu + dd)\P). The a term is 
renormalization point invariant in QCD, which implies that the matrix element in eqn. ( p8D 
will have an anomalous dimension of the same magnitude as that of the quark mass 
operator. Thus the zeroth moment of e{x) will change quite slowly with Q^, whereas 
we expect that the higher moments should scale approximately as 1/ \/Q^. Thus the 
distribution e(a;, Q^) will tend to increase at low x and decrease at high x as increases 

At twist-4 level there is no mixing of twist-2 and twist-3 operators with the local 



twist-4 operators [^, which simplifies the analysis somewhat. Of course there are a 
large number of local operators of twist-4 for each distribution function, and in the vast 
majority of cases their anomalous dimensions have not been calculated. Thus to evolve 
the twist-4 distributions f^{x) naive power counting, ~ is again used. In 

Figure 3 we see the twist-4 distributions at = 1 and 10 GeV^, where again the bag 
scale has been taken as /i = 0.4 GeV. 

Finally, the twist-two transversity distribution hi (x) will be of interest at scales above 



1 GeV . The anomalous dimensions of the twist-2 operators have been calculated so 
the evolution of the parton distribution from the bag scale is fairly straightforward. The 
results of this evolution up to = 1 GeV^ and = lo GeV^ are shown in Figure 4. 



^On completion of this work I learned of a calculation of the anomalous dimensions of e{x,Q^) 
which confirms these speculative remarks, in particular the anomalous dimensions of the operators cor- 
responding to the second and third moments of e are similar in magnitude to the anomalous dimensions 
of twist-2 distributions. 
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4 Nucleon Structure Functions 

The parton distribution functions calculated above can be measured in various processes, 
most notably deep inelastic lepton-nucleon scattering (DIS) and Drell-Yan (DY) processes. 
In this section we consider how the parton distributions combine to give the nucleon 
structure functions which can be measured in DIS, and those structure functions involving 
higher twist distributions which have been measured are compared with the calculated 
structure functions. 

In the Bjorken limit — >■ oo, x = Q"^ /2{P.q) fixed) the nucleon tensor W^^ can 



be expressed in terms of 14 structure functions |Tl 



V 



V 2v 



(f' ) ' ' ' ' V 

^q^q""-^-^ + (p^'q" +p''q^)^^ 
V 2v 

+M{Slf + Sir)^-^ + M^SIp"" + Slp^)^-^, (29) 

where the last seven structure functions (F3, 01^2,4,5, ^1,2) are related to parity-violating 
processes involving the weak interaction. I have introduced the shorthand notation p^ = 
P^ - ^q". 

It is conventional to introduce longitudinal structure functions, describing the scatter- 
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ing when the exchanged vector boson is longitudinally polarised 



Fl{x) = F2{x){l + ^—) -2xF^ix) 



iM'^x^^ 



/ AM'^x^\ 

aiix) = a2(s) ( IH \ - 2xai{x). (30) 



According to the Callen-Gross relation ||3^ both Fl and vanish in the Bjorken limit, 
and, ignoring QCD radiative corrections and nucleon mass effects, both are twist-4 struc- 
ture functions. 

The fourteen structure functions can all be expressed in terms of parton distribution 
functions |TI], O, |2H|. In the following all contributions from parton distributions with 



two or more light-cone fractions (ie distributions involving two quark fields plus one or 
two gluon fields and distributions involving four quark fields) will be dropped and the 
discussion is limited to the distributions calculated above. We also keep terms linear in 
the quark masses. The electroweak current of the quarks coupling to vector bosons 

is taken to be 

j'io = mri9v+9a75)m (31) 

where the vector and pseudo-vector couplings, and ga take the values of the standard 
model. Also, because weak currents can change quark flavour, the mass of a quark in the 
initial state is denoted as rrii, and the mass of a quark in the final state as mj. 

For unpolarised scattering there are five measureable structure functions. Fi and F2 (or 
Fi) describe electromagnetic scattering, F3 is related to parity violating weak scattering, 
and F4 and F5 are related to the scattering of non-conserved currents. We have: 

^ q 
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M 



(32) 



M 



F,(x) 



j:i-mM'+\9a,\')f? 



[X 



X 



(33) 
(34) 



where (—1)'^ is +1 for quarks and -1 for antiquarks. Finally 

M 



F^{x) 
F^ix) 



2Mx 



\^9vq^{m^ - mf) + \9aq^{mi + m/)] e^(x) 



H [lfl'f<?P("^i - ^f) + \9aq^{mi + m/) e'^{x). 



(35) 
(36) 



Thus Fl, F4 and F5 are twist-four structure functions. 

For scattering from a longitudinally polarised nucleon we have five structure functions: 

(37) 



1 



Glix) = -j:iM' + \9aq\')9Kx) 

^ q 



1 



ai[x) 



aL{x) 



T.i-mM'+\9aq\')9ti 



2 . 
2Mx 

M 



X] 



J2i~in\9vq\' + \9aq\') \-2Mx'glix) + 2m,xhlix) 



«4(X) = —J2i-'^yi\9vq\'^ + \9aq\'^)miXhl{x) 



a5{x) 



Q q 

2Mx 



[X] 



(38) 
(39) 
(40) 
(41) 



The twist-2 distribution ai is the longitudinally polarised analogue of F3, and if it were 
measured the role of the axial anomaly in the interpretation of measurements of Gi would 



be much clearer pTH . The structure functions 0x^4^5 are all twist-4, with 04^5 related to non- 
conserved currents. 

For scattering from a transversely polarised nucleon there are four structure functions: 



Gt{x) = -Y.{\9vq? + \9aq?)9'T{^) 



(42) 
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2M ^ 



bi{x) = Ei-mM' + \9aq\')xgUx) (44) 



b2{x) = ^j:i-^yiM" + \9aq\')mhUx). (45) 

These four structure functions are twist-3. Gt = Gi + G2 is related to the transverse 
spin of the quarks and antiquarks, 61 is its parity-violating partner, while G^ and 62 are 
related to non-conserved currents. 

Many experiments over the years have determined Fi, F3 and Gi to good precision, and 
in the near future Gt should be measured to similar precision. Fl has also been measured 
to a lower degree of precision, thus it should be possible for us to extract estimates of the 
parton distributions qt and f^. Unfortunately there is little hope of determining any of 



the other parton distributions from DIS. From the expressions for ^4^5, eqs. (|35| , |36|) one 
might hope that e{x) would be measurable in heavy quark production. However because 
\9vq\ = \gaq\ the terms in m/ cancel 0. A similar cancellation occurs in eqn. (|i3| ) making 
a determination of hi{x) very difficult. However the possibility exists to measure /ii,l(x) 
and e{x) in polarised Drell-Yan processes which may be possible at RHIC. 

Previous calculations [|T^ ^ have shown good agreement between experimental data 
for Fi{x) and F^lx) and the model predictions. We are now in a position to compare 
the model prediction for Fl{x) with experimental data. Experimentally the ratio R = 
Fl/2xFi has been measured [Q], and there is good evidence for a twist-4 contribution 



32| , Q at medium and high x and < 10 GeV^. Before we can compare with the data 



we should also take into account effects from perturbative QCD which give a non-zero 



■^This cancellation does not occur when the exchange boson is a Z boson, so there may exist a possibility 
to measure F4 5 at HERA 
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Fl at leading and next-to-leading order, and effects from the target nucleon mass. In the 
simplest version of the bag model the twist four contribution to R for electromagnetic 
scattering is given by 

In Figure 5 we show our calculated R^'^^{x,Q'^) at = 2 and 5 GeV^. We also show 
the experimental values of R and various perturbative QCD calculations, including target 
mass effects, and the effect of adding the twist-4 contribution from the bag model at 
these scales. Because R^^\x,Q'^) effectively scales as fi{x,jj?')/Q^, we calculate a large 
correction to the perturbative calculations at only low values of , however this twist-4 
correction is in the right direction, and for low gives a good fit to the data at medium 
X values. 

If the data at low improves in accuracy, it may be possible to perform an indepen- 
dent check on our value of the bag scale, /z, from a fit of our calculated parton distribution 
/4 with the data at various values of Q^. 

The other higher twist structure function that has been measured is G2{x) = Gt{x) — 
G\{x) |]. In Figure 6 the experimental data for G^ix) is compared with the calculated 
structure function at = 5 GeV^. The calculation is in reasonable agreement with 
the data. In particular the calculated G2[x) crosses the x axis in the right region. As 
this crossing is determined by the interplay between the twist-3 part of G2{x)i given by 

(3) 1 

gj- (x) and the twist-2 portion —gi{x) + dygi{y)/y, this gives us some confidence in our 
understanding of the relations between the twist-2 and twist-3 distributions. It will be 
interesting to check our model calculations against experimental at other values of Q^. As 
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increases we predict that the crossing point xq where ^2(2:^0) = will move to larger 



X. 



It is also interesting to note that while the bag model calculation for the structure 
function Gi{x) is not particularly good, presumably because the calculation does not 
include effects arising from the axial anomaly [0, in the case of G2{x) the agreement 
with the data is better, probably because of the cancellation between the two terms in 
the twist-2 part of G2{x). 

Finally in Table 1 the calculated values for some of moments of the twist-2 and twist-3 
parts of G2{x) are displayed and compared with values from other calculations and those 
extracted from the experimental data. 

Before concluding, there are some limitations to these calculations which should be 
borne in mind. There is no flavour dependence of the calculated parton distributions, 
so the distributions and structure functions shown above should be taken as applying to 
scattering from an isoscalar target. For our comparison of R this is not important, as 
experiment has shown no target dependence for R but for G2 the flavour dependence 
will be as important as for Gi[13, 37|. In the model, the flavour dependence of the parton 



distributions comes from the SU (6) spin-flavour part of the wavefunction, and the one 
gluon exchange mechanism means that the two-quark intermediate state must be either 
a scalar (S=0) or vector (S=l) with masses and > m^, respectively |3^, |13 



Furthermore no account has been taken of contributions to the parton distributions from 
intermediate states consisting of four quarks. However these contributions have only a 



small effect on valence distributions, and this is all in the small x region |T3[ , which is not 
of great importance here. 
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An important correction to the parton distributions will come from mesonic contribu- 



tions via the so-called Sullivan process |3^, where meson {tt, p, K . . .) exchange between 



the virtual photon and the nucleon occurs. These processes are known to lead to impor- 
tant modifications in the parton distributions at twist-2 particularly for gi{x), 
and there is no reason to suppose that the higher twist distributions should be unaffected. 
This should be important in the model calculation of the structure function G2- 

5 Summary 

In this paper we have discussed the various parton distribution functions up to twist four 
in the context of the MIT bag model. The discussion has been limited to those distribu- 
tions involving two-quark correlation functions, however these include effects from the bag 
boundary, which is a simple model of the effects of the gluon field. From the twist-3 distri- 
butions, which involve quark-gluon correlations, we have calculated the twist-3 structure 
function, G2{x), and found the calculation in good agreement with the experimental data, 
which gives us some confidence in using the bag as a model of the quark-gluon dynamics 
involved in DIS. 

Four-quark correlation functions are expected to be small as they involve the cal- 
culation of Hill- Wheeler-type overlap functions for four quarks in the bag. Also these 
distributions are not directly accessible to experiment. However in future work we shall 
look at the corrections four-quark distributions make to the two-quark distributions cal- 
culated here. 

At twist-four we have seen that the calculated distribution f4{x) gives a contribution 
to the experimentally measured R = Fl/2xF2. When added to the contributions from 
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perturbative QCD and target mass effects, the calculation gives good agreement with the 
data at (5^ ~ 1 GeV^ an dmedium values of x. 

Some improvements can be made to our calculations. On the model side we can add 
flavour dependence to our distributions, and we can also add the important effects of the 
pion cloud. Also we can improve our evolution techniques for the distribution functions. 
For the twist-2 distributions next-to-leading-order QCD corrections are available. At 
twist-3 there have been recent developments in the calculations of anomalous dimensions 
of the relevant operators, so that a correct leading order QCD approach is possible, though 
this will be very difficult. 

I would like to thank the Institute for Theoretical Physics, University of Adelaide, for 
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Figure Captions 

Figure 1. The parton distribution functions calculated at the bag scale = fj,"^. 
(a) The twist-two distributions fi{x),gi{x) and hi{x). (b) The twist-three distributions 
e{x),hL{x) and gxix). (c) The twist-four distributions f4,{x),g3{x) and h^{x). 
Figure 2. The twist-three parton distributions at (a) = 1 GeV^ and (b) = 10 
GeV2. 

Figure 3. The twist-four parton distributions at (a) = 1 GeV^ and (b) = 10 
GeV2. 

Figure 4. The transversity dependent twist-two parton distribution hi{x,Q'^) at = 1 
GeV2 and = iq GeV^. 

Figure 5(a). The twist-four distribution R^^\x^Q'^) calculated at = 2 GeV^ and 
Q2 = 5 GeV^. 

Figure 5(b). Experimental data on i? and QCD calculations with (upper band) and 
without (lower band) the twist-four contribution R^^\x) at = 2 GeV^. 
Figure 6. Comparison of experimental data |^ on G^{x) with this calculation. 
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M,[G^,] 




This work 


-5.9 X 10-3 


-1.2 X 10-3 


-0.3 X 10-3 


CM. Bag Model 38| 


-1.2 X 10-3 






Lattice QCD|| 


(-26.1 ±3.8) X 10-3 






Data§ 


(-6.3 ± 1.8) X 10-3 


(-2.3 ±0.6) X 10-3 


(-1.0 ±0.3) X 10-3 







< 


"6 


This work 


6.3 X 10-3 


1.3 X 10-3 


0.3 X 10-3 


CM. Bag Model 


-17.4 X 10-3 






Lattice QCD 


(-48 ± 5) X 10-3 






QCD Sum Rule 4C] 


(-6 ±3) X 10-3 






Datag 


(5.4 ±5) X 10-3 


(0.7 ± 1.7) X 10-3 


(0.1 ±0.8) X 10-3 



Table 1. Comparison of various calculations with experimental data for the moments of 
Gl and the twist-3 matrix elements < = 2(^A^„[G'f] ± Mn[G{]) = 2^A<„[Gg^^]. 
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Fig. 1(a) Twist-2 Parton Distributions 




Fig. 1(b) Twist-3 Parton Distributions 




Fig. 1(c) Twist-4 Parton Distributions 
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Fig. 2(a) Twist-3 Parton Distributions 
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Fig. 2(b) Twist-3 Parton Distributions 
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Fig. 3(a) Twist-4 Parton Distributions 




Fig. 3(b) Twist-4 Parton Distributions 



0.05 



0.04 



0.03 - 



0.02 - 



0.01 












0.2 



0.4 



0.6 



0.8 



f_4(Q^2=10) 
g_3(Q^2=10) 
h_3(Q^2=10) 



X 



Fig. 4 Twist-2 Trans versity Distribution 




Fig. 5(a) Twist-4 Contribution to R 
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